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Abstract. This paper considers extreme values attained by a centered, multidimen- 
sional Gaussian process X(t) = (-Xi(f), . . . , X n (t)) minus drift d(t) = (di(t), . . . , d n (t)), 
on an arbitrary set T. Under mild regularity conditions, we establish the asymptotics of 



logP ( 3t e T : P| {Xi(t) - di(t) > q>u} 

V i=i 



for positive thresholds qi > 0, i = 1, . . . ,n and u — ► oo. Our findings generalize and 
extend previously known results for the single-dimensional and two-dimensional cases. A 
number of examples illustrate the theory. 



1. Introduction 

Owing to its relevance in various application domains, in the theory of stochastic processes, 
substantial attention has been paid to estimating the tail distribution of the maximum 
value attained. In mathematical terms, the setting considered involves an M-valued sto- 
chastic process X = {X(t) : t £ T} for some arbitrary set T and a threshold level u > 0, 
where the focus is on characterizing the probability 



(1.1) Pj^supX(t) >uj =¥(3teT : X(t) > u) . 

More specifically, the case in which X is a Gaussian process has been studied in detail. This 
hardly led to any explicit results for (jl.ip . but there is quite a large body of literature on 
results for the asymptotic regime in which u grows large. The prototype case dealt with a 
centered Gaussian process with bounded trajectories for which the logarithmic asymptotics 
were found: it was shown that 

(1.2) lim u~ 2 logP ( supX(t) > u ) = - (2<4) \ 
where 

o\ := sup EX 2 (t). 

teT 

See Adler jj], p. 42] or Lifshits 0, Section 12] for this and related results. The monographs 



Lifshits [9( and Piterbarg llj contain more refined results: under appropriate conditions, 
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an explicit function 4>(u) is given such that the ratio of (jl.ip and 4>(u) tends to 1 as u — > oo 
(so-called exact asymptotics) . The logarithmic asymptotics (|1.2|) can easily be extended 
to the case of noncentered Gausssian processes if the mean function is bounded. The 
situation gets interesting if both trajectories and the mean function of the process are 
unbounded. In this respect we mention Dufneld and O'Connell 0] and Debicki jsj, where 
the logarithmic asymptotics of P(sup t>0 (X(i) — d(t)) > u) for general centered Gaussian 
processes X, under some regularity assumptions on the drift function d, were derived; see 
also Hiisler and Piterbarg Dieker 0] and references therein. 

While the above results all relate to one-dimensional suprema, considerably less attention 
has been paid to their multidimensional counterparts. One of few exceptions is provided 



by the work of Piterbarg and Stamatovic i2j], who considered the case of two R- valued 



possibly dependent, centered Gaussian processes {Xi(t\) : t\ £ Ti} and {X 2 (t 2 ) : t 2 £ T 2 }. 
They found the logarithmic asymptotics of 

(1.3) P(3(ti, t 2 ) € T : Xtih) > u, X 2 (t 2 ) > u) 

for some T C T\ x T 2 , under the assumption that the trajectories of X\ and X 2 are 
bounded. 

In this paper our objective is to obtain the logarithmic asymptotics of (following the 
convention that vectors are written in bold) 

(1.4) P(u) :=w(3t€T: f]{X t (t) - d t (t) > q^ ; 

here {X(t) : t € T}, with X(t) = (Xi(t), . . . ,X n (t))', is an ]R n -valued centered Gaussian 
processes defined on an arbitrary set T C R m , for some m, n S N, the d{{-) are drift 
functions and qt > are threshold levels, i = 1, . . . ,n. Our setup is rich enough to cover 
both of the cases in which P(u) corresponds to the event in which (i) it is required that 
there is a single time epoch t E M such that Xi(t) — di(t) > q^u for all i = 1, . . . ,n and 
(ii) there are n epochs (ii, . . . , t n ) such that Xi{ti) — di(ti) > for all i = 1, . . . , n. We 
get back to this issue in detail in Remark [lj where it is also noted that the theory covers 
a variety of situations between these two extreme situations. 

Compared to the one-dimensional setting, the multidimensional case requires various tech- 
nical complications to be settled. The derivations of logarithmic asymptotics usually rely 
on an upper and lower bound, where the latter is based on the inequality 

P(u) > supP (f]{Xi(t) - ck(t) > q iU }) . 

^ T V'=l / 

Strikingly, in terms of the logarithmic asymptotics, this lower bound is actually tight, 
which is essentially due to the common 'large deviations heuristic': the decay rate of the 
probability of a union of events coincides with the decay rate of the most likely event 
among these events. A first contribution of the present paper is that we show that this 
argument essentially carries over to the multidimensional setting. In order to obtain the 
lower bound one needs asymptotics of tail probabilities that correspond to multivariate 
normal distributions. In this domain a wealth of results are available (see, e.g., Hashorva [?J 
and references therein), but for our purposes we need estimates which are, in some specific 
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sense, uniform. A version of such estimates, that is tailored to our needs, is presented in 
Lemma 01 

The upper bound is based on what we call a 'saddle point equality' presented in Lemma 
[TJ It essentially allows us to approximate suprema of multidimensional Gaussian process 
X by a specific one-dimensional Gaussian process, namely a properly weighted sum of 
the coordinates Xi of X. Formally, we identify weights Wi = Wi(t,u) > such that the 
inequality 

(n n \ 

3teT: J2mMt) > ^2m{uqi + ck(t)) , 
i=l i=\ / 

is logarithmically asymptotically exact, as u — > oo. The reduction of the dimension of 
the problem allows us to use one-dimensional techniques (such as the celebrated Borell 
inequality). Interestingly, the optimal weights can be interpreted in terms of the solution 
to a convex programming problem that corresponds to an associated Legendre transform 
of the covariance matrix of X. A different weighting technique has been developed in 
Piterbarg and Stamatovic 12j for the case n = 2, but without a motivation for the weights 
chosen. We recover the result from 12| in Remark [5j Our analysis of (jl.4p extends the 
results from 0, El], in the first place because M"- valued Gaussian processes are covered 
(for arbitrary n £ N). The other main improvement relates to the considerable generality 
in terms of the drift functions allowed; these were not covered in |l2j |. 
The paper is organized as follows. In Section [2] we introduce notation, describe in detail 
objects of main interest to us, and state our main result; we also pay special attention to 
the rationale behind the assumptions that we impose. In Section [3] we illustrate the main 
theorem by presenting a number of examples; one of these relates to Gaussian processes 
with regularly varying variance functions. We also explain the potential application of our 
result in queueing and insurance theory. In Section|3]we describe how the multidimensional 
process X can be approximated by a one-dimensional process Z, obtained by appropriately 
weighting the coordinates X{. We prove some preliminary results about the characteristics 
of the process Z . This section also contains the saddle point equality mentioned above, 
Lemma [JJ which is the crucial element of the proof of our main result. Section 0] also 
contains all other lemmas needed to prove Theorem [TJ as well as the proof of our main 
result itself. 



2. Model, notation, and the main theorem 

In this section we formally introduce the model, state the main theorem, and provide the 
intuition behind the assumptions imposed. 

2.1. Model and notation. Let T C M m , for some m G N. In this paper we con- 
sider an W 1 - valued (separable) centered Gaussian process X = {X(t),t £ T} given 
by X(t) = (Xi(t), . . . , X n (t))'. Let the so-called drift function be denoted by d(t) = 
(di(t), . . . ,d n (t))'. Now, denote the covariance matrix of X(t) by S^. Throughout the 
paper it is assumed that the matrix T,t is invertible for every t G T. Here and in the 
sequel, we use the following notation and conventions: 

• We say v > w if V{ > W{ for all % = 1, . . . , n. 

■ We write diag(i>) for the diagonal matrix with V{ on the ith position of the diagonal. 

• We define vw := diag(v)w' = (v±wi, . . . , v n w n )'. 
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• For a G R, we let i(a) be an n-dimensional vector (a,..., a)' and also let = 

(o,...,oy. 

• We adopt the usual definitions of norms of vectors := ((x, x)) 1 / 2 , where (-, •) 
is the Euclidean inner product. 

• We let f(u) ~ g(u) denote that lim u _>oo f(,u)/g(u) = 1. 

• We write R™ := {x G R n : x > 0,x 0}. 

Throughout the paper not all vectors are of dimension n (for instance t is of dimension 
m), but the above notation should be understood with obvious changes. 
With each we associate the matrix Kt = (kij(t))ij< n , defined as 

K t = dia g (a 1 ; 1 1/2 (t) ; . . . , a-v 2 (t))E t - 1 di ag (9 1 ; 1 1/2 (t), . . . , a$ 2 (t)) 

with T,^ 1 = (dij(t))ij< n . We mention that kij(t) G [—1, 1] and that —kij(t) is commonly 
interpreted as some sort of partial correlation between Xi(t) and Xj{t) controlling all other 
variables X/%(i), k ^ i,j. 

2.2. Main result. Throughout the paper, we impose the following assumptions. 
Al sup tGT hj(t) < 1 for all £ ^ j, i, j = 1, . . . , n. 

A2 sup ter (Xj(£) — edi(t)) < oo a.s. for all i = 1, . . . , n and all e G (0, 1]. 

If a process X and a drift function d comply with assumptions A1-A2, then to shorten 
the notation, we will write that (X,d) satisfies A1-A2. 
For a point t G T and a vector q > 0, define 

M XAq (u,t) := inf (v + d^^iv + dit))), 

v>uq 

M XAq (u;T) := - inf M XAq (u,t). 

With these preliminaries we are ready to state our main result. The following theorem can 
be seen as an n-dimensional extension of [12, Theorem 1] and [5,, Theorem 2.1]. 

Theorem 1. Assume that (X,d) satisfies A1-A2. Then, for any q > 0, 

(2.1) logP(3t G T : X(t) - d(t) > uq) M XAq (u;T) as u -> oo. 

Remark 1. The result stated in Theorem [1] enables us to analyze, with Tj C R, 

(2.2) P f f] { sup (Xi(ti) - di(U)) > m 

V'=l Ueii 

To see this, let T := T\ x . . . x T n . Also define processes \Yi(t) : t G T}, i = 1, . . . , n, such 
that Yi(t) := X^(tj), for i = 1, ... ,n. Analogously, let mj(t) := di(ti), i = 1, ... ,n. Then 
(1231) equals 

P(3t G r : Y(t) - m(t) > ng) , 

which, under the proviso that Al- A2 are complied with by the newly constructed (Y, m), 
fits in the framework of Theorem [1] This example naturally extends to the situation where 
the sets Tj are of dimension higher than 1. 
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2.3. Discussion of the assumptions. In this subsection we motivate the assumptions 
that we imposed. 

Remark 2. Assumption Al plays a crucial role in the proof of Lemma [H It can be 
geometrically interpreted as follows. For a fixed t G T, the distribution of X(t) equals 
that of Bt M, where Bt is a matrix such that T,t = BfB' t and M is an M n -valued standard 
normal random variable. For some quadrant Qt, we need in the proof of Lemma 2] a lower 
estimate of ¥(X(t) G Q t ) = F(M G B^Qt). For i = 1, ...,n let be, as usual, the 

j 



standard basis vectors of R n . Then the cosine of the angle ati j between B t e.{ and B t 1 e 



is given by 

C0S(Qij) _ <V-^>> = (^') = . _ Mt) . 

*" ||BT 'cllllB, - 'ejll K <=illlW r'ejll v'*; (')««(<) 

We thus observe that Al entails that, for all i G T, there is no pair of vector B^ l ei and 
B^ 1 ei, with j 7^ j, that 'essentially coincide', i.e., the angles remain bounded away from 
0. Therefore, for any a; G B^ 1 Qt, one can always find a set At such that a; G At C B^ l Qt 
and At has a diameter that is bounded, and a volume that is bounded away from zero, 
uniformly in t G T. 

Remark 3. For e = 1, assumption A2 assures that the event 

\J{X(t)-d(t)>uq} 

teT 

is not satisfied trivially. The following example shows that if A2 is not complied with, 
then it is not ensured that we remain in the realm of exponential decay. Consider a one- 
dimensional case in which X = {X (t) : t > 0} is a standard Brownian motion, and for any 
5 > let d(t) := (1 + S)y/2t log logt. From the law of the iterated logarithm we conclude 
that the process X does not satisfy A2 for every e G (0, 1]. On the other hand we have 
(take t := u A ) 



(sup (x(t) - (1 + 5) ^2t log logt) > u J > P f 

Vt>o v ' / yi + (i 



> 1 



+ 5)uy / 21og(41ogi 

where here is the real-valued standard normal random variable. On the logarithmic 
scale the latter probability behaves roughly, for u large, as 

-(1 + 5) 2 log logu. 

For the case of n = 1, A2 has been required in 0, Theorem 2.1] as well. 

Remark 4. The drift functions dj, i = l,...,n, are not assumed to be increasing, but 
under assumption A2 we have ti := initeT di(t) > — oo. Because we are interested in the 
asymptotic behavior of the probability in (|2.ip as u — > oo, we can assume that u > uq := 
— mmi(£i / qi) , and therefore the coordinates of uq + d(t) stay positive for all t G T. In 
what follows we shall always assume that u > Uq. 
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3. Examples 

In this section we present examples that demonstrate the consequences of Theorem [TJ We 
focus on computing the decay rate Mx : d,q( u ]T) in two cases: (i) the case of X having 
bounded sample paths a.s.; (ii) the case of the Xi having stationary increments, regularly 
varying variance functions, and di(-) being linear. While in the former example the drift 
functions do not influence the asymptotics, in the latter example the drifts do have an 
impact on the decay rate. 

3.1. Bounded sample paths and drift function. We here analyze the case of (X,d) 
satisfying 

Bl The process X has bounded sample paths a.s. 

B2 There exists D < oo such that < D for all t G T and i = 1, . . . , n. 

We note that under B1-B2, it trivially holds that assumption A2 is complied with as well. 
Assumptions B1-B2 are satisfied when T is compact, X has continuous sample paths a.s. 
and d is continuous for instance. Let us introduce the following notation 

I X , q {T) := inf inf (v^v) . 

' y teTv>q x 1 ' 

The following corollary is an immediate consequence of Theorem [TJ 

Proposition 1. Assume that (X,d) satisfies Al and B1-B2. Then, 

u 2 

logP(3t G T : X(t) - d(t) > uq) —I x , q {T), as u -> oo. 

The above proposition states that in the 'bounded case' that we are currently considering, 
we encounter the same asymptotic decay as in the driftless case (d = 0). 

Remark 5. Some special cases of Proposition [TJ have been treated before in the literature. 
In particular, let X\ = {X\(t\) : t\ £ T±} and X2 = {^2(^2) : ^2 S ^2} be two centered and 
bounded M- valued Gaussian processes. We introduce the notation <7j(tj) := y / Var(A"j(tj)), 
r(t) := Corr(Ai(ti),X 2 (t 2 )) and also 

qi a 2 (t 2 ) (Ti(ti) q 2 



c q (t) :-- 



mm 



.<ri{h) 12 Qi o- 2 (t 2 ) 
Then, upon combining Proposition [JJ with Remark [JJ we obtain, with T C T\ x T 2 , 

logP(3(t 1 ,t 2 ) G T : Ax (t0 > q lU ,X 2 (t 2 ) > q 2 u) 

„ -± inf \ ( 1+ Mt)-r(t))> 

2 (*i,ta)er {m.in{o- 1 (t 1 )/q 1 ,o- 2 {t 2 )/q 2 }) 2 \ 1 - r 2 (t) 

as u — y 00. Observe that the above formula is also valid for r(i) = ±1. This recovers the 
result of Piterbarg and Stamatovic [13] • 

3.2. Stationary increments, linear drift. This section focuses on the logarithmic asymp- 
totics of {X(t) — i(t) : t > 0}, where X(t) = SY(t) for some invertible matrix S and, as 
usual, Y{t) = (Yi(t), . . . ,Y n (t)y. We assume that, for i = 1, . . . , n, 

CI {Yi(t) : t > 0} are mutually independent, M-valued, centered Gaussian processes 
with stationary increments. 
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C2 The variance functions a 2 (t) := Var(Yj(i)) are regularly varying at oo with indexes 
G (0, 2). Without loss of generality we assume that < or < . . . < a„ < 2. Moreover, 
assume that there exists k £ {1, ... ,n} such that of ~ . . . ~ c K a 2 for some q > and 
lim^oo a K (t) I a K+ i(t) = (if k = 1, then set c K = 1; if the first condition is satisfied with 
k = n, then the second one is redundant). 

C3 lim t ^o &i {t)\ log |t|| 1+£ < oo for some e > 0. 
We analyze 

(3.1) P (3i > : X{t) - i(t) > uq) . 

Probabilities of this type play an important role in risk theory, describing the probability 
of simultaneous ruin of multiple (dependent) companies; see Avram et al. [3] for related 
results. The one-dimensional counterpart of (|3.1|) was considered in Debicki 0] in the 
context of Gaussian fluid models. Related examples and further references can be found 
in the monograph 10]. In the following proposition we derive the logarithmic asymptotics 

of (prm . 

With q as in C2, set 

C := diag(l,c 2 , . . . , c K , 0, . . . , 0) 

and 

J(C, S, q, a) := mf mf . 

Proposition 2. Assume that Y satisfies C1-C3, and S is an invertible matrix. Then, 
for {X(t) :t>0}:= {SY(t) :t>0}, 

u 2 

logP(3t > : X(t) - i(t) > uq) r j(C,S,q,ai), as u -> oo. 

2a f (u) 

Proof. We start by checking that A1-A2 are satisfied for (X,i). Indeed, let us note that 
the matrix Kt = K is constant. Besides, since S is invertible, then K is invertible too, 
which combined with the fact that K is positive-definite and fc^, = 1, straightforwardly 
implies that assumption Al is satisfied. 

Since Y has stationary increments, then under C1-C3 lim^oo Yi(t)/t = almost surely 
and therefore (using that X consists of linear combinations of the Yi, i = 1, . . . ,n) as- 
sumption A2 is complied with; see 0, Lemma 3] for details. Now following Theorem 

m 

Mx,i,,(«;[0,oo)) = iinf inf + i(t)), ^S'^v + *(*))> 

Z t>0 v>uq 

= \ inf inf (S-\uv + ui{t)),R-JS-\uv + ui(t))) 

L t>0 v>q 

= % inf inf <5- 1 (^ + i(t)), J R- i 1 S- 1 (^ + i(t))>, 

Z t>0 v>q 

where the matrix R^ 1 equals diag(<7j~ 2 (i), . . . , a~ 2 (t)), which is the inverse of the covariance 
matrix of Y. Using the regular variation of of (■), we find that, as u — > oo, 

a 2 (u)R~t -> t~ ai C, as u -)■ oo. 
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By virtue of the uniform convergence theorem we arrive at 

m t \h w u2 ■ f f (s-Hv + i(t)),cs-Hv + i(t))) 

M x ^u- [0,oo)) ~ ^mf mf , 

as u — > oo. This completes the proof. □ 

4. The proof of the main theorem 

This section is devoted to the proof of our main result - Theorem [TJ We will achieve this 
by establishing an upper bound and a lower bound. We start by presenting the following 
'saddle point equality' that plays a crucial role in the upper bound. 

Lemma 1. Let A be any positive- definite matrix. Then, 



(w,q) 2 . , i . - 
sup / , = mf (v,A l v) 

{W,AW) v>q" 



for any vector q £ M". Moreover, if v* is the optimizer of the infimum problem in the 
right-hand side, then w* := A~ 1 v* is an optimizer of the supremum problem in the left- 
hand side. 

Proof. Decompose A = BB' for some nondegenerate matrix B. Then, 

(w,q) 2 (w,q) 2 , x . x 2 

— and Iv.A v) = \\B v . 



(w,Aw) \\B'w\\ 2 
Now, for w £ Wt, the Cauchy-Schwarz inequality yields 

(w,q) = inf (w,v) = inf (B'w,B~ 1 v) < \\B'w\\ inf ||£> _1 t;||. 

Dividing both sides by ||-B'iw|| > and optimizing the left-hand side of the previous display, 
we arrive at 

sup < w 'g> 2 < inf (v^A-'v). 

{W,AW) v>q 

To show the opposite inequality, assume that v* is such that 

inf (v^A^v) = (v*,A- l v*). 

v>q 

The Lagrangian function of the above problem is given by L(v, A) := (v, A~ 1 v} — (A, v — q) 
for A > 0, and due to complementary-slackness considerations we necessarily have that 
A~ l v* > 0, and if (A^v^i > 0, then v*i = Thus take w* = A~ l v* £ W% , so that 

(„•„>» (A-^, q f v 
(w*,Aw*) {A- l v*,v*) x 1 
Indeed, the last equality is equivalent to 

(A- 1 v*,q-v*) = Q, 
but recall that if (A^v*)^ ^ 0, then (q — v*)i = 0. Hence finally, 

su p r^>™ i ( v ' A ~ lv )> 

w( z R n (W,AW) v>q 

which proves the opposite inequality. This finishes the proof. □ 
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The main idea behind the proof of the upper bound of Theorem [1] is that the M n -valued 
process X(t) — d(t) can be effectively replaced by a suitably chosen M-valued Gaussian 
process. The asymptotics of the latter process can then be handled using the familiar 
techniques for real- valued Gaussian processes. 
For any vector w G WL, define 



Zu,w \t) 



(w, uq + d[t)) ' 
and observe that (with u > uq; cf. Remark H|) 

P(3t G T : X{t) - d(t) > uq) < P ( sup Z uw {t) > 1 J . 

\teT ' / 

The vector w in the process Z U)W can be seen as a vector of weights assigned to the 
coordinates of X. For fixed u and w the process Z u ^ w is a centered Gaussian process. We 
shall show that it also has almost surely bounded sample paths. 

Lemma 2. Under A1-A2, the process Z uw is a centered Gaussian process with bounded 
sample paths almost surely, for each w £ Wl and u > uq. Moreover, 

p 

sup Z UjW (t) — > as u — > oo. 
teT 

Proof. Without loss of generality we can assume that ||u;|| = 1. For any L > 1, recalling 
the definition of I from Remark [H 



supZ nw (i) >L) =¥(3teT: (w,X(t)) > (w, Luq + LI + L(d(t) - £))) 

< P (3t G T : (w,X(t)) > {w, L(uq + £) + (d(t) - £))) 

< P (3t G T : (u>, X(t) - d(t)) > (to, L(uq + £)-£)) 

< P l^2wiSup(Xi(t) - ck(t)) > L(w,uq + £) - (w,£)\ 



vi=l teT 



< 



\^sup(Xj(t) - di(t)) + > Lmm(uqi + li)j\fn 
^ teT * 



vi=l 



where the last probability tends to zero with L — >• oo due to A2. This proves that Z UyW 
has bounded sample paths almost surely. 

The last probability also tends to zero with L > 1 fixed and u — >• oo. On the other hand, 
for any L < 1 we have 

supZ uu; (t) > L ) = P (3t G T : (tO,X(t) - Ld(t)) > L{w,uq)) 
teT ' / 



< P Vsup(Xi(t) - Ldi(t)) + > uL mingi/yn , 



vi=l * eT 



where the last probability tends to zero with u oo by virtue of A2. We therefore have 
that sup tgT ^ u ,io (t) converges to in probability. □ 
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The above considerations remain true even if w depends on u and t. This observation 
allows us to optimize the variance of the process Z UjW , while retaining its sample path 
properties. Notice that 

(w,uq + d(t)) 

Therefore, take w* = w*(u,t) such that 

(4 1) (w*,S t n>*) _ . nf (w,-E t w) 

(w*,uq + d(t)) 2 {w,uq + d(t)) 2 

and denote by Y u (t) the process Z u ^ w *(t) with the weights w = w* chosen as above. Let 
o- 2 (t) be the variance function of the process Y u (t). Then, by Lemma [fl 

(4.2) a- 2 (t) = M x ,d, q (u,t). 

To estimate the tail of the supremum of the process Y u (t) we intend to use Borell's in- 
equality [l], Theorem 2.1]. To apply this result, we need to verify that the expectation of 
su PtgT Y u {t) vanishes as u — > oo. This is done in the next lemma. 

Lemma 3. Under A1-A2, with uq as in Remark^ 

(1) Mx,d.q(u;T) > for each u > uq; 

(2) lim^oo M x ,d, q (u;T) = oo; 

(3) lim u ^ 00 Esup tgT y u (t) = 0. 

Proof. From Lemma [2] we know that for a fixed u the process Y u has bounded sample 
paths almost surely. This implies that sup tgr °u(t) < °°- But 

supa 2 u (t) = suviMx.a^t))- 1 = Umx^u-T))- 1 
teT teT z 

and claim (1) follows. 

The proof of (2) is a consequence of the fact that under A2 



sup Y u (t) > 1 J — > as u — > oo, 
teT J 

and for M being a standard normal random variable 



su P y u (t) > 1 > supP (Y u (t) > 1) = P [AT > inf JM XAq (u,t) . 
teT J teT V teT v / 

To prove the last claim, observe that the almost sure boundedness of sample paths of Y u (t) 
implies that E sup teT Y u (t) < oo and it easily follows that the family (sup tgT Y u (t)) u is 
uniformly integrable. Now claim (3) follows from the second part of Lemma [2j □ 

Before we proceed to the proof of Theorem [1] we state a technical lemma, which is a 
prerequisite for the proof of the lower bound. 

Lemma 4. Under Al, there exist constants C\ < oo, C2 > such that for any t € T 
logP(X(i) - d(t) > uq) > -\M XAq {u,t) - CiMi /2 dq («,t) + C 2 . 
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Proof. Set 

Qt := {x G R n : x > uq + d(t)}, 

and let B t be such that B t B' t = E t . Then X(t) = B t Af, where N is an K n -valued standard 
normal random variable with the density function 

f{x) = £> n exp ( -- (x,x) 



2 

for some normalizing constant D n . In this notation, we have 

P (X(t) - d(t) >uq)=F (X(t) G Q t ) = P (M G B^Qt) • 

Now let x* = x*(u,t) G B^ l Q t be such that 

M x ,d, q {u,t) = inf (aj.E^a;) = inf (as, as) = (sc*,sc*) , 
x ^Qt xeB 1 r l Qt 

and let A t := B(x*, 1) n B~ l Q t , where B(x*, 1) is a ball in M n of radius 1 and center x* . 
Then, 

P (N G B^Q*) > / /(su) dx. 
Set A (a;, 3:*) := {x, x) - (x*, x*) . Then 

P [M G B^Qt) > B n Vol(^ t )exp f- ]-M XAq {u, t) - \ sup A(x,x*)) . 

\ 1 1 xeA t J 

Since 

A(x, x*) < 2\\x - x*\\ (x*, x*) 1/2 + - x*\\ 2 , 

we have that 

sup A(35,35*) < 2diam(A t )M^ /2 d (u,t) + diam 2 (^ t ). 
xeA t ' ' 

Therefore the claim follows if diam(At) and Vol(-At) can be bounded uniformly in t G T 

from above and below, respectively. 

Observe that, by the construction of At, diam(At) < 1. Besides, the quadrant Qt is 
spanned by the standard basis (ej) in 1" fixed in the point uq + d(t). The cosine of the 
angle atij between B^ 1 ej and B t r 1 e J - is given by cos(ajj) = fcjj; see Remark[2j Under Al 
this angle is bounded away from zero, uniformly in t G T. Therefore inft^T Vol (.At) > 0. 
This completes the proof. □ 

Now we are ready to prove the main theorem. 

Proof of TheoremUl Put P(u) := P (3t G T : X(t) - d(t) > uq). We split the proof into 
two parts: the lower and the upper bound. 

Lower bound: The lower bound follows directly from Lemma [Hand the inequality 

logP(u) > suplogP(X(i) - d(t) > uq) . 
teT 

Upper bound: Let w* : M + x T — > Wt be the mapping chosen in (|4.ip . Now as in the 
definition of the process Y u , 

P(u) < P (3t G T : (to*, X(t)) > (to*, uq + d(t))) 

(w*,X(t)) \ ( 

sup \ 2m*\\ > 1 = M su P y « * > 1 
teT (w*,uq + d(t)) J \ t£T 
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where the passage from the re-dimensional quadrant to the tangent increases the probabil- 
ity. Recall that the variance a* (t) of Y u (t) equals (M XA? (M,t))^; cf. (02]). Moreover, 
thanks to Lemma O the Gaussian process Y u has bounded sample paths almost surely. 
Therefore, Borell's inequality implies that 

supY u (i) > l] < 2exp f - (l-EsupY u {t)] M x , d , q (u;T) ] . 
vteT / \ \ t£T / J 

Now from (2) and (3) of Lemma [3] we obtain 

,. iogP(su Pter r u (t) > l) 

hmsup — — ^- — — < —1 

M ^oo M x ,d,q{u;T) 

and the claim follows. □ 

Remark 6. From the proof of the upper bound we obtain the useful inequality 

P (3t G T : u?*X(t) > io*(ug + d(t))) < P (3t € T : (to*, X(t)) > {w*, uq + d(t))) , 

which we have proven to be exact in terms of logarithmic asymptotics. Let v* = v*(u,t) 
be such that 

(v* + d(t),Zt\v* + d(t))) = inf (v + d(t),Zt\v + d(t))) . 

v>uq 

Then the optimal weights w* are given by w*(u,t) = T l ^ t 1 v*(u,t), or alternatively, due to 
Lemma [H by 

(w,uq + d(t)) 2 

w (u,t) = arg sup — r . 

{w, T. t w) 

Observe that the weights do not depend on u in the case of d = 0. 
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